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Using an approximate time-dependent density functional theory method, we cal- 

S: 

■ culate the absorption and luminescence spectra for hydrogen passivated silicon 

nanoscale structures with large aspect ratio. The effect of electron confinement 
in axial and radial directions is systematically investigated. Excited state relaxation 
leads to significant Stokes shifts for short nanorods with lengths less than 2 nm, but 
has little effect on the luminescence intensity. The formation of self-trapped excitons 
is likewise observed for short nanostructures only; longer wires exhibit fully delocal- 

c\i ; 

ized excitons with neglible geometrical distortion at the excited state minimum. 
1 PACS numbers: 
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> : I. INTRODUCTION 

Research on silicon nanowires (SiNW) has intensified in the past years, owing to their po- 
tential applications in future nano-technologies, such as nanosensors^, nanoscale electronics 
and photonics devices.-^ This interest results from the electronic structure of SiNW beeing 
critically dependent on the size, orientation, passivation and doping level of the nanostruc- 
ture. For example, it has been shown experimentally that the band gap can be tuned by 
choosing different growth directions and diameters for the wire.- On the theoretical side, the 
ground state electronic structure and transport properties of different kinds of SiNW have 
been investigated by several authors .-^^i^i^ 



2 



Light absorption and emission of these systems have also attracted great attention, be- 
cause the band gap in bulk crystalline Si is small and indirect, while that in SiNW can 
become large and direct due to the quantum confinement effect.— iMiUy&±2 This experimen- 
tal progress paves the way to obtain visible light from silicon materials, and provides for 
the possibility of future optoelectronics applications. In the past, theoretical descriptions 
of the optical properties of SiNW were mainly restricted to determining their absorption 
spectrum.-ii&ii To quantify luminescence, it is often tacitly assumed that de-excitation oc- 
curs resonantly. However, absorption and emission spectra can differ considerably due to 
excited state relaxation, and simulations that take this Stokes shift into account are currently 
missing for SiNW. 



II. METHOD AND SIMULATION APPROACH 



In this letter, we use the time-dependent density-functional based tight-binding (TD- 
DFTB) method^ (for a recent review see^), to study the excited states and optical prop- 
erties of SiNW. The same method has been applied to silicon quantum dots,— i2ii22i2iL2ii2£i2& 
where it was found to provide a high degree of reliability and computational efficiency com- 
pared with the parental time-dependent density functional theory (TD-DFT) from which it 
is derived. For example, the optical gap of Si5Hi 2 is predicted to be 6.40 eV, in very good 
agreement with the experimental value of 6.50 eV.— Also for larger clusters like Si35H36, the 
TD-DFTB estimate for the lowest allowed singlet transition (4.37 eV) coincides with high 
level ab-initio results from multi- reference second order perturbation theory (4.33 eV).— 

In principle, SiNW should be treated as quasi one- dimensional periodic systems, since 
their lengths are generally in the micrometer range while their diameters are several nanome- 
ters only.™ However, the TD-DFTB method we use here is currently restricted to deal with 
finite systems, thus, we construct silicon nanostructures along the (110) direction with in- 
creasing length to approach the experimentally realized SiNW. This allows us to demonstrate 
how the localized excitons in silicon quantum dots evolve into delocalized ones as the length 
of the Si nanostructures continously increases. Focussing on finite structures also guarantees 
that we do not leave the trust region of TD-DFT. For periodic systems, local or semi-local 
approximations to the exchange-correlation functional lead to a collapse of the many-body 
excited state energies on the ground state Kohn-Sham gap.— ^ A better description of the 



3 




FIG. 1: (color online) Ball and stick models of four (110) silicon nanowires (SiNW) with different 
diameters viewed from the top. Here Si atoms are yellow-colored and H atoms are grey-colored. 

electron-hole interaction is provided by quasi-particle calculations at the GW level combined 
with the solution of the Bethe-Salpeter equation.— 1 ^ For the confined systems discussed here, 
however, TD-DFT provides a reliable and affordable way to investigate optical properties. 

We studied (110) SiNW with four different diameters, as shown in[TJ The diameters d 
are estimated to be 0.84 nm, 0.98 nm, 1.08 nm, 1.24 nm, respectively. Due to the limitation 
of computational resources, SiNW with other growth directions or larger diameters were not 
considered. The simulation protocol for each structure can be summarized as follows: First, 
to build the structural model, we fully saturated the dangling bonds of surface Si atoms with 
hydrogen atoms. Second, the constructed models were relaxed in the ground state by con- 
jugate gradient optimization using the DFTB method. 33 Third, we obtained the low energy 
part of the absorption spectrum at the ground state optimum by solving for the energies 
and oscillator strengths for the lowest five singlet excitations. And fourth, assuming rapid 
internal conversion, (i.e. the validity of Kasha's rule), the emission energies are evaluated 
by geometry optimization in the first excited singlet state (Si). For most systems studied 
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here, this assumption does not represent an approximation, since the strongest absorbing 
state is actually the Si. In the optimization process, we took advantage of a recent imple- 
mentation of analytical excited state forces for the TD-DFTB scheme.— In all calculations 
the basis set consists of s and p orbitals for Si atoms and an s orbital for H atoms, and the 
gradient-corrected PBE exchange-correlation functional is employed 35 . 

For a given structure, excited state energies (ujj) are obtained in two steps.— ^ First, 
ground state spin-restricted DFTB calculations are performed to obtain the Kohn-Sham 
(KS) orbitals ipi and the KS energies 6j. These single-particle values are corrected following 
the TD-DFT linear response treatment of Casida 36 

(1) 

klr 

where o^- = £j — and a and r are spin indices. The coupling matrix Kij a>k i T defined as 3 ^ 

K ^ kW = II ^ (r) ^ (r)( F^ + ^(r)^(rO ) 
xip k (r')ipi{r')drdr' ', 

is further simplified in the TD-DFTB approach using the Mulliken approximation.—^ 
With the results of equation (CQ), oscillator strengths are calculated as 3 ^ 

/' = h E i£< > \/^(4t + 4i)i 2 ( 2 ) 

■j , — V u i 

k=x,y,z ij ' 

III. RESULTS AND DISCUSSION 

In El we show the absorption and luminescence spectra for the thinnest (110) nanowire 
(d = 0.84 nm) with different lengths, from I = 2.1nm (Siu2H 98 ) to I = A.2nm (Six^Hi^) . 
Similar to what we have previously observed in silicon quantum dots with increasing 
diameter,— <2i we find here that both absorption and emission energies slightly red-shift 
with increasing wire length. Moreover, a sizable Stokes shift of around 0.1 ~ 0.2 eV is 
discernible. Since excited state relaxation does not lead to a significant change in oscillator 
strengths, all considered nanostructures should exhibit significant luminescence intensity. 

Next, we extracted the energies of the first allowed transition for SiNW with different 
diameters to investigate the size dependence of the optical gap (see [3]). When the length 
of the Si nanostructures increases from 0.39 nm to 2.73 nm, the absorption energies 
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FIG. 2: (color online) Absorption (solid lines) and emission spectra (dashed lines) for (110) 
nanowires with diameter 0.84 nm and different lengths (from 2.7 nm to 4.3 nm). The spectra 
have been broadened by 0.01 eV to simulate finite temperature. 

will decrease monotonously from around 5 eV to around 3 eV. This is due to the decrease 
of quantum confinement effects. Concurrently, the emission energies E em i show an overall 
increase (although with some fluctuations for most cases) from 1 eV to around 3 eV, and 
then follow the trend of the absorption energies. This general behaviour has also been found 
in similar studies on the size dependence of the optical gap in spherical Si quantum dots.— >^ 
Quantum confinement effects are also observed in the radial dimension. Considering a fixed 
length of I = 1.56 nm, the optical gap decreases from 3.66 eV for the thinnest diameter 
(d = 0.84 nm) to 3.21 eV for d = 1.24 nm. 

In order to understand the transition from quantum dots to wires more deeply, we analyze 
the geometrical distortions arising from the excited state relaxation. In small quantum dots, 
one Si-Si bond is extremely stretched up to 2.70 A.— Depending on the actual structure, the 
location of this bond can be either in the center of the cluster or at its surface. When the rods 
are long enough, the geometry distortions are small and more homogenously distributed over 
all bonds. For the thinnest nanostructures {d = 0.84 nm), the distortions locate at one Si-Si 
bond in the center of the structures for lengths less than 1.95 nm and delocalize beyond this 
value, owing to the confinement provided by the rigid surrounding layer.— For a diameter of 
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FIG. 3: TD-DFTB absorption energies (circles) and emission energies (squares) for (110) SiNW 
with different diameter and lengths. 

1.08 nm, geometry changes are observed for one Si-Si bond at the surface of the structures 
up to a length of 1.17 nm (Si^Hqs). The situation for the other two cases (d = 0.98 nm 
and d = 1.24 nm) is more involved. For the wire with d = 0.98 nm, distortions are found 
primarily for a Si-Si bond at the center of the quantum dots below 0.78 nm (Si^H^s). The 
distortion is then delocalized for the cluster Si^H^ at 1.17 nm, which corresponds to the 
peak in the emission energy curve in [3j It should be mentioned here that our simulation 
protocol locates the local minimum closest to the the Frank-Condon point, i.e. the initial 
structure after light absorption. It cannot be ruled out that the global minimum of the Si 
potential energy surface is again of localized nature. Whether luminescence occurs also from 
this state will crucially depend on the excited state lifetime and the kinetic energy of the 
ions after absorption. Excited state molecular dynamics simulations could shed light on this 
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interesting question, which is however outside the scope of this study. 

Turning back to the results, a surface Si-Si bond is found to be elongated for the cluster 
SigoHgo at 1.56 nm and thereafter all distortions have a delocalized character for longer 
structures. A similar case are the wires with d = 1.24 nm, where a delocalized distortion 
appears for SIqqHqq (I = 0.78 nm), and a localized distortion for SigoH 76 (/ = 1.16 nm). 

Additional information on the electronic structure can be extracted from the KS molecular 
orbitals that are involved in the light absorption and emission.— We take the rods Si^H^ 
{I = 1.17nm) and Sig^I!^ (I = 2.32 nm) with diameter d = 0.84 as examples. In both 
systems, the Si excited state wavefunction is dominated by a single-particle transitions 
from the highest occupied molecular orbital (HOMO) to the lowest unoccupied molecular 
orbital (LUMO). In HI these orbitals are depicted for the optimized geometries in So and Si . 
We find for the ground state minimum of both structures that the HOMO distributes mainly 
along the central zone of the wire, while the LUMO is located on the surface. The molecular 
orbitals at the excited state minimum conformation are quite different for the two structures 
under investigation. The shorter structure, Si4$H 5Q , exhibits a LUMO concentrated at the 
center of the cluster that gives rise to a repulsive force. The significant Si-Si bond stretch in 
combination with the resulting sizable Stokes shift is completely in line with the self-trapped 
exciton model of Allan, Delerue, and Lannoo.— ^ In contrast, the LUMO of the longer rod, 
SiQ 6 H 8 Q, distributes more or less homogenously along the structure, indicating a delocalized 
exciton and geometry distortion. 



IV. CONCLUDING REMARKS 



In conclusion, we simulated the absorption and emission spectra of SiNW with different 
diameters. The evolution of the optical properties from small quantum dots to nanowires 
with large aspect ratio has been investigated. While short nanorods with lengths below 
2 nm show localized excitations and the formation of self-trapped excitons, excited state 
relaxation has little effect on longer structures which exhibit delocalized excited states. 
Notwithstanding these general trends, silicon nanostructures with similar extensions can 
exhibit quite different localization characteristics conjoined with largely differing emission 
profiles. This fact may partially explain the reported spread of photoluminescence energies^ 
and calls for atomistic simulations of these systems which take this strong conformation 
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FIG. 4: (color online) Absolute values of the molecular orbitals HOMO (blue) and LUMO (red) for 
ground state optimized (left) and excited state optimized (right) structures of Si^H^q (up) and 
SiqqHsq (down), which are models of (110) SiNW with diameter d = 0.84 nm. The plot corresponds 
to an isovalue of 0.001. 

dependence into account. 
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